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Abstract 

We have measured quantum transport through an individual Fe4 single-molecule magnet em- 
bedded in a three-terminal device geometry. The characteristic zero-field splittings of adjacent 
charge states and their magnetic field evolution are observed in inelastic tunneling spectroscopy. 
We demonstrate that the molecule retains its magnetic properties, and moreover, that the magnetic 
anisotropy is significantly enhanced by reversible electron addition / subtraction controlled with 
the gate voltage. Single-molecule magnetism can thus be electrically controlled. 
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Rationally designed magnetic molecules [TJ [2] can be used as building blocks in future 
nanoelectronic devices for molecular spintronics [3J, classical [4] and quantum information 
processing (51 E] • They usually have long spin coherence and spin relaxation times due to 
a weak spin-orbit and hyperfme [7] coupling to the environment. Of crucial importance for 
applications is the ability to adjust the magnetic properties by external stimuli. In bulk 
samples tuning magnetic properties by light has already been demonstrated [8], but on the 
single-molecule level it has not been achieved. In addition, to tune magnetic properties on a 
single-molecule level, the use of local electric fields is preferred as it allows for a direct and 
fast spin-state control. 

Addressing individual magnetic molecules on chip [9] has proven to be extremely chal- 
lenging. Attempts to incorporate archetypal single-molecule magnets (SMMs) Mni2 into 
a three-terminal device [TO] ITT] have been followed by observations that these complexes 
undergo electronic alterations when self-assembled on gold surfaces [12]. In this letter we 
demonstrate electric-field control over the magnetic properties of an individual Fe4 molecule 
connected in a planar three-terminal junction. In the neutral state the bulk properties of 
this SMM are well documented [13J [T4] and more importantly they are retained upon depo- 
sition on gold [15]. From transport measurements we find that a Fe4 molecule in a junction 
can still behave as a nanoscale magnet with the anisotropy barrier close to the bulk value. 
In addition, upon reduction or oxidation induced by the gate voltage the barrier height 
increases; i.e., by charging, the molecule becomes a better magnet. 

Characteristic of a SMM is its magnetic anisotropy, creating an energy barrier U which 
opposes spin reversal; i.e., the high spin of the molecule points along a preferred easy axis, 
making it a nanoscale magnet. The anisotropy lifts the degeneracy of the ground high- 
spin multiplet, even in the absence of a magnetic field. The splitting of the lowest two 
levels is known as the zero-field splitting (ZFS); see Figure lb. The magnetic anisotropy 
is described by the parameter D, which for the Fe4 molecule in the bulk phase equals 
D = 0.051—0.056 meV; the ground state spin for this class of molecules S — 5, the anisotropy 
barrier U = DS 2 £ 1.3 - 1.4 meV and the ZFS is (25 - 1)D = 0.46 - 0.50 meV [13] [14]. 

We use Fe4 molecules with formula [Fe4L2(dpm)g] (Hdpm = 2,2,6,6-tetramethyl-heptan- 
3,5-dione) (see Figure 1). Two derivatives, Fe 4 Ph and Fe 4 CgSAc, were synthesized 
by functionalizing the ligand H 3 L = R-C(CH 2 OH)3 with R = phenyl and R = 9- 
(acetylsulfanyl)nonyl, respectively, and were prepared as described elsewhere [T31[T4]. Their 
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stability in a dry toluene solution at a 8 mM concentration was checked using spectroscopic 
techniques (see Supporting Information). Nanometer-spaced electrodes were fabricated us- 
ing self-breaking [TB] of an electromigrated gold wire in a toluene solution of the molecules at 
room temperature [17J. Three-terminal transport measurements were performed at a tem- 
perature T = 1.6 K. By varying the gate voltage V g the molecular levels are shifted, thereby 
providing access to magnetic properties in adjacent charge states of the SMM. Currently 
this is not possible with other techniques and therefore little is known about the magnetic 
properties of charged SMMs. 

To quantify the anisotropy of an individual molecule we have performed transport- 
spectroscopy measurements in the Coulomb blockade regime. In total we have measured 
648 devices, of which 48 showed Coulomb blockade signatures and two-dimensional con- 
ductance maps pjZ] (dl/dV versus gate and bias voltage) were measured on these. The 
observation of molecule-related features (e.g. excitations) depends strongly on the domi- 
nant transport mechanism and electronic coupling to the leads. Typically, the strength of 
the electronic coupling in molecular junctions is significant so that the resolution of single- 
electron tunneling spectroscopy is limited by the tunnel-broadening. In this case sharper 
higher-order inelastic cotunneling peaks [18] may resolve the ZFS, but if their broadening 
due to Kondo correlations exceeds the ZFS, only a single, broad zero-bias Kondo peak re- 
mains, masking the low-bias ZFS features. We identified 9 stable devices with indications 
of transport through a high-spin molecule: four devices showed broad Kondo peaks in two 
adjacent charge states, indicating that at least one of the charge states is a high-spin state. 
These Kondo features will be the subject of a subsequent paper. Two samples showed a 
magnetic field-dependent transition at an energy scale of 1.5 meV and another one showed 
a low-bias current suppression indicative of spin blockade [19], involving a high-spin ground 
state. The microscopic origins of these observations are not clear at the moment; it should 
be noted, however, that the 1.5 meV value is close to the anisotropy barrier U. Two devices 
showed transitions in the inelastic cotunneling spectra at an energy scale below 1 meV, close 
to the bulk value of the ZFS. In this paper we focus on these two devices. 

In Figure 2a and 3a,b we present the measurements as differential conductance maps 
at zero magnetic field. Sample A (Figure 2) features the Fe4Ph derivative and sample B 
(Figure 3) the Fe^gSAc derivative. The plots show conducting regions characteristic of 
single-electron tunneling (SET). Outside the SET-regime, the molecule is in a well-defined 
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charge state with N electrons and the current is suppressed (Coulomb blockade). Figure 2a 
and 3a,b clearly show a sizable conductance in the blockade regime due to tunnel processes 
of second ( "cotunneling" ) or higher order. Here, inelastic cotunneling dl/dV steps or peaks 
appear at a bias voltage V& = ±A/|e|, and their observation allows for the determination of 
molecular excitation energies A in a particular charge state [18]. The charge state can be 
changed by the application of a gate voltage when crossing a highly conductive SET region. 
By making the gate voltage more positive, the electron number increases and the molecule 
is reduced; by making it more negative, the electron number decreases and the molecule is 
oxidized. 

In Figure 2a (sample A) three inelastic cotunneling lines are visible in the left charge 
state. For large negative gate voltages these occur at V = ±0.6, ±4.6 and ±6.7 mV. For the 
right charge state there are lines at T4 — ±0.9 mV and ±5 mV. For sample B cotunneling 
excitations appear at V = ±0.9 mV and ±5 mV (Figure 3a) in the left charge state and at 
V = ±0.6 mV and ±7 mV (Figure 3b) in the right one. 

Additional information on the nature of the low-energy excitations can be obtained by 
measuring their evolution in a magnetic field. In Figure 2b, c and Figure 3c, d dl/dV is 
plotted as a function of V for different magnetic field values at two adjacent charge states 
for each sample. The energy of the lowest excitation increases with magnetic field without 
splitting and is symmetric upon field reversal (see Supporting Information). Such behavior 
is a hallmark of ZFS described by the spin Hamiltonian [21 12"UH2"2"] 

H = -DS 2 z + 9f i B B-~S, (1) 

where the first term is the uniaxial magnetic anisotropy with the easy axis z. The second 
term is the Zeeman interaction of spin S with magnetic field B, where g is the Lande factor 
and \xb the Bohr magneton. Importantly, model (1) predicts a non-linear dependence of the 
excitation energy on the magnetic field if the angle between the easy axis of the molecule 
and the magnetic field is substantial (see Supporting Information). This non-linearity is 
caused by the mixing of \S Z ) states due to the presence of a transverse component of the 
field. 

To quantitatively compare the data with the model, the excitation energy has been 
determined from individual dl/dV curves. For sample A (Figure 2b, c), we have taken 
the peak positions as the ZFS excitation energy and this energy is plotted versus magnetic 
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field in Figure 2d,e. For sample B, the inelastic cotunneling excitations have the form of 
conductance steps, reflecting a weaker molecule-electrode tunneling T than for sample A. 
We have fitted the measured dl/dV to a Lambe-Jaklevic formula [231 IM] and plotted the 
excitation energies in Figure 3e,f. 

We first compare the data of sample A with the model. For the left charge state (Fig- 
ure 2b) the energy of the first excitation at B = T is close to the ZFS value (0.5 meV) in 
the bulk [13] . We assume that this is the neutral state, and since the Fe4 molecule maintains 
its ground state spin S = 5 when deposited on gold [IS], we estimate D = 0.06 meV, so 
that U = 1.5 meV. The increase of the excitation energy with B is linear, implying an 
angle 9 between B and the easy- axis below 45°. It should be noted that the angle cannot 
be controlled in our junctions and is expected to vary for different samples. Using model 
(1) we have calculated the energy difference between the ground and the first excited state 
originating from the same side of the anisotropy barrier as a function of applied magnetic 
field. Using the angle 9 as an additional parameter, the best fit of the data to the model, 
shown in Figure 2, yields g = 2.1 and 9 = 0°. 

Reduction and oxidation of the molecule inevitably change its magnetic properties, al- 
though it is not a priori clear in which way. Three-terminal spectroscopic measurements 
can provide an answer since the ZFS and the change in the total spin upon charging can 
be obtained independently. First, the difference in spin values of adjacent charge states can 
be determined from the shift of the degeneracy-point in a magnetic field (see Supporting 
Information). We infer an increased S = 11/2 in the reduced charge state of sample A. 
From the measured ZFS we then find an enhanced D = 0.09 meV and U = 2.7 meV. From 
the magnetic field dependence, we infer g = 1.8 and 9 = 0° (Figure 2e) when fitting the 
data to the model (1). For sample B we observe the bulk ZFS value in the right charge state 
(Figure 3d,f) and therefore identify this as the neutral state with S = 5 and D = 0.06 meV. 
A clear non-linear Zeeman effect is now observed, indicating that the field is at a substan- 
tial angle with the easy-axis. We have fitted the data to the model obtaining g — 2.1 and 
9 = 71°. Also for the left charge state the Zeeman effect is non-linear (Figure 3e) and a 
three-parameter fit yields D = 0.09 meV, g = 2.0 and 9 = 60°. Here, we have used S = 9/2 
for the oxidized state (see Supporting Information) and estimate U = 1.8 meV. Note that for 
small magnetic fields the data show a deviation from the model, which is not yet understood. 

In view of the rich SMM excitation spectrum the question arises why only the ZFS ex- 
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citation is observed up to a bias of several mV. We have performed extensive calculations 
which indicate that this indeed should be the case (see also Supporting Information). In Fig- 
ure 4a,b we show dl/dV maps calculated using quantum kinetic equations (KE), accounting 
for tunnel processes to first and second order in T [25]. We model the low-energy spectrum 
of two successive charge states including the charge-dependent ZFS. Figure 4 shows that for 
the experimental temperature the ZFS dl /dV-step is dominant, showing a linear Zeeman 
effect for 9 = 0° (Figure 4a) and a non-linear behavior for 9 = 71° (Figure 4b). Numerical 
renormalization group (NRG) calculations using a Kondo model [26J cover the regime where 
the tunnel coupling dominates over the thermal energy (in contrast to the KE approach). 
Figure 4d,e show that the ZFS now appears as a peak due to the exchange scattering through 
the SMM indicating significant Kondo correlations. 

In summary, we have demonstrated electric-field control over the anisotropy of a single 
magnetic molecule in a three-terminal junction. We found a stronger magnetic anisotropy 
upon both reduction and oxidation induced by the gate voltage. This enhancement may be 
related to the alteration of single-ion anisotropy, which should be substantial when changing 
the redox state, but more studies including quantum chemistry calculations have to be per- 
formed to confirm this. Our findings open the route to new approaches in tuning magnetism 
on a molecular scale and the possibility of manipulating individual magnetic molecules for 
future use in nanoelectronic applications. 

This work was supported by NanoSci-ERA, FOM, DFG SPP-1243. 

Supplementary information containing chemical characterization data, 
transport measurements and calculations details can be found at 



http://pubs.acs.org/doi/suppl/10.1021/nll009603 free of charge. 



[1] Kahn, O. Molecular Magnetism; VCH Publishers, Inc.: New York, 1993. 
[2] Gatteschi, D.; Sessoli, R.; Villain, J. Molecular Nanomagnets; Oxford University Press: Ox- 
ford, 2006. 

[3] Bogani, L.; Wernsdorfer, W. Nature Materials 2008, 7, 179-186. 

[4] Affronte, M. J. Mater. Chem. 2009, 19, 1731-1737. 

[5] Leuenberger, M. N.; Loss, D. Nature 2001, 410, 789-793. 



6 



[6] Lehmann, J.; Gaita-Arino, A.; Coronado, E.; Loss, D. Nature Nanotech. 2007, 2, 312-317. 
[7] Bertaina, S.; Gambarelli, S.; Mitra, T.; Tsukerblat, B.; Muller, A.; Barbara, B. Nature 2008, 
453, 203-206. 

[8] Herrera, J. M.; Marvaud, V.; Verdaguer, M.; Marrot, J.; Kalisz, M.; Mathoniere, C. Angew. 

Chem. Int. Ed. 2004, 43, 5468-5471. 
[9] Grose, J. E.; Tarn, E. S.; Timm, C; Scheloske, M.; Ulgut, B.; Parks, J. J.; Abrufia, H. D.; 

Harneit, W.; Ralph, D. C. Nature Materials 2008, 7, 884-889. 
[10] Heersche, H. B.; de Groot, Z.; Folk, J. A.; van der Zant, H. S. J.; Romeike, C; Wegewijs, M. R.; 

Zobbi, L.; Barreca, D.; Tondello, E.; Cornia, A. Phys. Rev. Lett. 2006, 96, 206801. 
[11] Jo, M.-H.; Grose, J. E.; Baheti, K.; Deshmukh, M. M.; Sokol, J. J.; Rumberger, E. M.; 

Hendrickson, D. N.; Long, J. R.; Park, H.; Ralph, D. C. Nano Lett. 2006, 6, 2014-2020. 
[12] Mannini, M.; Sainctavit, P.; Sessoli, R.; Cartier dit Moulin, C; Pineider, F.; Arrio, M.-A.; 

Cornia, A.; Gatteschi, D. Chem. Eur. J. 2008, 14, 7530-7535. 
[13] Accorsi, S.; Barra, A.-L.; Caneschi, A.; Chastanet, G.; Cornia, A.; Fabretti, A. C; Gat- 
teschi, D.; Mortalo, C; Olivieri, E.; Parenti, F.; Rosa, P.; Sessoli, R.; Sorace, L.; Wernsdor- 

fer, W.; Zobbi, L. J. Am. Chem. Soc. 2006, 128, 4742-4755. 
[14] Gregoli, L.; Danieli, C; Barra, A.-L.; Neugebauer, P.; Pellegrino, G.; Poneti, G.; Sessoli, R.; 

Cornia, A. Chem. Eur. J. 2009, 15, 6456-6467. 
[15] Mannini, M.; Pineider, F.; Sainctavit, P.; Danieli, C; Otero, E.; Sciancalepore, C; Ta- 

larico, A. M.; Arrio, M.-A.; Cornia, A.; Gatteschi, D.; Sessoli, R. Nature Materials 2009, 

8, 194-197. 

[16] O'Neill, K.; Osorio, E. A.; van der Zant, H. S. J. Appl. Phys. Lett. 2007, 90, 133109. 

[17] Osorio, E. A.; Bj0rnholm, T.; Lehn, J.-M.; Ruben, M.; van der Zant, H. S. J. J. Phys.: 

Condens. Matter 2008, 20, 374121. 
[18] De Franceschi, S.; Sasaki, S.; Elzerman, J. M.; van der Wiel, W. G.; Tarucha, S.; Kouwen- 

hoven, L. P. Phys. Rev. Lett. 2001, 86, 878-881. 
[19] Weinmann, D.; Hausler, W.; Kramer, B. Phys. Rev. Lett. 1995, 74, 984-987. 
[20] Sessoli, R.; Gatteschi, D.; Caneschi, A.; Novak, M. A. Nature 1993, 365, 141-143. 
[21] Hirjibehedin, C. F.; Lin, C.-Y.; Otte, A. F.; Ternes, M.; Lutz, C. P.; Jones, B. A.; Hein- 

rich, A. J. Science 2007, 317, 1199-1203. 
[22] Tsukahara, N.; Noto, K.; Ohara, M.; Shiraki, S.; Takagi, N.; Takata, Y.; Miyawaki, J.; 

7 



Taguchi, M.; Chainani, A.; Shin, S.; Kawai, M. Phys. Rev. Lett. 2009, 102, 167203. 
[23] Lambe, J.; Jaklevic, R. C. Phys. Rev. 1968, 165, 821-832. 

[24] Kogan, A.; Amasha, S.; Goldhaber-Gordon, D.; Granger, G.; Kastner, M. A.; Shtrikman, H. 

Phys. Rev. Lett. 2004, 93, 166602. 
[25] Leijnse, M.; Wegewijs, M. R. Phys. Rev. B 2008, 78, 235424. 

[26] Romeike, C; Wegewijs, M. R.; Hofstetter, W.; Schoeller, H. Phys. Rev. Lett. 2006, 96, 196601. 
[27] Carretta, S.; Santini, P.; Amoretti, G.; Guidi, T.; Caciuffo, R.; Candini, A.; Cornia, A.; 
Gatteschi, D.; Plazanet, M.; Stride, J. A. Phys. Rev. B 2004, 70, 214403. 



8 



a 




FIG. 1: Fe4 molecule, (a) Structure of the Fe4 molecules (color code: iron=purple, oxygen=red, 
carbon=grey, sulfur =yellow) . Left: Magnetic core with four S = 5/2 iron (III) ions antiferromag- 
netically coupled to give a molecular spin S = 5. Center: Fe^gSAc derivative. Right: Fe4Ph 
derivative, (b) Energy diagram of the ground spin multiplet at zero field. The S z ^ levels 
corresponding to different orientations of the spin vector along the easy axis of the molecule are 
doubly-degenerate. The S z = +5 and S z = —5 states are separated by a parabolic anisotropy 
barrier U. An important property of the Fe4 molecule is the large exchange gap to the next S = 4 
high-spin multiplet [27] in the neutral state, which is 4.80 meV and 4.65 meV, for Fe4Ph and 
Fe4CgSAc, respectively |13|. 114]. Transport below a bias voltage of a few mV therefore only probes 
magnetic excitations of the ground high-spin multiplet, in contrast to the Mni2 derivatives |10j . 
(c) Drawing of a three-terminal junction with a single Fe4Ph molecule bridging two gold electrodes 
(yellow) on top of an oxidized aluminum gate (grey). 
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FIG. 2: Characteristics of sample A. (a) Color plot of dl/dV versus V and Vq at T = 1.6 K and 
B = T. In the black regions in the middle of the plot single-electron transport (SET) is allowed. 
The lock-in amplifier saturates in these high conductance regions, (b) dl/dV as a function of V for 
Vg = -1-5 V and various magnetic field values. Successive curves are offset by 250 nS. (c) Same as 
(b) for Vq = 2 V and an offset of 400 nS. (d) and (e) Excitation energy as a function of magnetic 
field for the same Vq values as in (b) and (c). Red lines are fits with D = 0.06 meV, g = 2.1, S = 
5 in (d) and D = 0.09 meV, g = 1.8, S = 11/2 in (e) and 6 = 0° for both cases. 
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FIG. 3: Characteristics of sample B. (a) and (b) Color plots of dl/dV versus V and Vq at T = 
1.6 K and B = 0T highlighting the behavior in two adjacent charge states, (c) dl/dV as a function 
of V at Vq = V for successive magnetic fields B. Curves are offset by 125 nS. Red lines are fits 
to a Lambe-Jaklevic formula (see Supporting Information), (d) Same as (c) for Vq = 1.6 V with 
curve offsets 50, 80, 140 and 220 nS. (e) Energy of the first excitation as a function of B at Vq = 
0.1 V. The red line is a fit to model (1) for S = 9/2, D = 0.09 meV, g = 2.0 and 9 = 60°. (f) Same 
as (e) for V G = 1.1 V. The red line is a fit for S = 5, D = 0.06 meV, g = 2.1 and 9 = 71°. 



11 



dl/dV (a.u.) b dl/dV (a.u.) 




FIG. 4: Calculated dl/dV as a function of bias voltage and magnetuc field, (a) KE result for 8=0° 
using the parameters estimated for sample A with electron temperature T = 1.6 K and the gold 
conduction bandwidth W = 8.1 eV. Shown is the field evolution of the inelastic cotunneling step 
for the 5 = 5 state at gate voltage Vq = -1.5 V. The conductance is scaled to its maximum value, 
(b) Same as (a) but using the parameters estimated for sample B with gate voltage Vq = 1-6 V 
and 9 = 71°. (c) dl/dV traces taken from (b) corresponding to sample B. (d) NRG result for 8=0° 
for the Kondo model using the same parameters as in (a) with zero temperature and an exchange 
tunneling constant J = 0.1/ p (p is density of states), (e) Same as (d), but now for 9 = 71°. (f) 
dl/dV traces taken from (d) corresponding to sample A. 
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